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1. Introduction 

Hercules X-l (discovered by Tananbaum et al. 1972) is one of the most frequently 
observed X-ray binary systems. The intermediate mass of the donor star, 2.2 M , leads 
to a wealth of behavior seen in both low mass and high mass systems. Although we suspect 
the mass flow is primarily through Roche lobe overflow giving rise to an accretion disk, as in 
the Low Mass X-ray Binaries (LMXB), the variable P Cygni profiles observed in UV lines 
suggest either a transient stellar wind or a stellar wind that is photoionized in some regions 
(Boroson, Kallman, & Vrtilek 2001). There are theoretical grounds as well to expect that 
winds may arise in this system when X-rays from the neutron star heat the surface of the 
normal star or accretion disk (Arons 1973; Davidson & Ostriker 1973; Basko et al. 1977; 
London, McCray, & Auer 1981; Begelman, McKee, & Shields 1983; Begelman & McKee 
1983). 

Most LMXBs, including Z-sources and atoll-sources, do not show persistent pulsations, 
perhaps because they have neutron stars with low magnetic fields. These sources thus lack an 



-2 - 



essential window into their kinematics. Light travel-time delays in the 1.24 second Her X-l 
X-ray pulsation period determine that the neutron star is in a nearly circular 1.7 day orbit 
with semimajor axis a x = 13.86 seconds. The orbital period is slowly lengthening from mass 
loss in the system (Deeter et al. 1991). 

Fortuitously, Her X-l has high orbital inclination (> 80 degrees), allowing total X-ray 
eclipses to constrain the size of the donor star, and the high galactic latitude offers only small 
reddening, allowing the system to be observed at crucial UV and even EUV wavelengths. 

The system's most mysterious variability is the 35-day X-ray high and low cycle (Giac- 
coni et al. 1973). For 8-11 days (the "Main-On state") of this cycle Her X-l emits X-rays 
with a ~ 10 37 erg s -1 total X-ray luminosity. For a ~ 4 day cycle (the "Short-On state"), 
halfway through the 35 day cycle, Her X-l's output is several times lower. For the remaining 
portion of the 35 days, the Off state, the observed X-ray flux is only several percent of that 
during the Main-On state. 

The behavior of the optical flux and X-ray spectrum over the 35-day period shows that 
the 35-day variability is not isotropic. Instead, the X-ray emission is merely obscured by a 
warped accretion disk whose shape "precesses" globally throughout a 35-day period. The 
optical emission from the normal star varies with the 1.7 day binary orbit. Its spectral type 
is variously classified as A through F, because X-rays from the neutron star heat the surface 
of the facing side of the Roche lobe. The spectral type continues to change with the orbital 
period throughout the longer 35-day cycle, implying that X-rays continue to heat the surface 
of HZ Her. 

The transition between the Off and Main-On state is only a few hours, whereas the 
Short-On state is entered into more gradually. This indicates that regions of the disk at 
different radii may obscure the central source, in which case the disk must precess globally 
with the same period. The variation in pulse shape with the long-term period is sometimes 
taken to imply that the innermost edge of the accretion disk precesses with the 35 day period, 
obscuring local regions on the neutron star (Scott, Leahy, & Wilson, 2000). 

The physical cause of a 35-day global disk precessing disk warp, or a definitive relation 
between system parameters and disk precession period, remains unknown and is a major goal 
of investigations of the system. It has been proposed, but not confirmed, that reprocessed 
X-ray radiation pressure (Maloney & Begelman 1997; Wijers & Pringle 1999) or an X-ray 
driven disk wind and corona (Meyer, Meyer-Hofmeister, 1984) could propagate the warped 
shape in the disk. Similar "long-term periods" have been observed in other X-ray binaries, 
including LMC X-4. Her X-l serves as a prototype of this behavior. 

X-ray dips present another variability that has not been explained definitively. Crosa 
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& Boynton (1980) showed that the dips recur every 1.65 days, near, but not at, the 1.62 day 
beat period between orbital and precessional periods. This was confirmed by long-term X- 
ray observations (Scott & Leahy, 1999). The dips are thought to be associated with the gas 
stream between the stars, but there is no consensus explanation for them or their period. 

As befits such a frequently observed, enigmatic, yet prototypical source, the methods of 
spectroscopy have been applied to the system at a wide range of wavelengths. These studies 
have borne fruit with determinations of the system's elemental abundances and the orbital 
motion of the normal star (the neutron star's pulse delays indicate only the motion of one 
component of the binary). 

High resolution spectroscopy in the X-ray range using XMM (Jimenez-Garate et al. 
2002) shows a multitude of X-ray emission lines presumably from an accretion disk corona. 
The line ratios indicate that the gas is enhanced from CNO processing from a massive 
progenitor. Jimenez-Garate et al. (2005) confirmed the CNO enhancement by observations of 
more than two dozen emission lines originating in an accretion disk corona. Model predictions 
of the disk corona's response to illumination by the central X-ray source are in reasonable 
agreement with the observed fluxes for low and moderate Z elements ( O through S), but the 
Fe xxv, Fe xxvi lines are several times brighter than predicted (Jimenez-Garate et al. 2005). 

From observations of optical absorption lines (Reynolds et al. 1997), and from optical 
pulsations that result from the X-rays from the neutron star periodically striking the surface 
of the normal star (Middleditch & Nelson 1976), we know that the normal star has a mass 
~ 2.2 M and the neutron star has a mass 1.5±0.3 M . The optical spectrum shows 
absorption lines from HZ Her, emission at the Bowen blend near 4640A (see Schachter, 
Filippenko, & Kahn 1989 for a discussion of the formation of these lines), and emission from 
HellA4686. Still et al. (1997) performed Doppler tomography on these optical emission lines 
but could not conclude that they arose in a symmetric accretion disk. 

The UV spectrum is crucial for an understanding of the accretion disk and the X-ray 
illuminated face of HZ Her. The continuum emission from both the illuminated star and the 
disk peaks in the UV and the disk contributes a greater fraction than at optical wavelengths 
(Cheng, Vrtilek, & Raymond 1995). From models of the variable UV continuum as observed 
with the International Ultraviolet Explorer (IUE, Vrtilek & Cheng 1996) showed that a 
change in the accretion disk precession could explain an anomalously low period of X-ray 
emission. 

To study global accretion one should observe gas at ionization stages and temperatures 
resulting from X-ray illumination of the disk, and one would need to resolve Doppler-shifted 
velocities that correspond to the orbital motion of the neutron star (160 km s -1 ), motion in 
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the accretion disk (expected to be p=j 300 km s 1 at the edge of the disk), and the velocities 
expected in a disk or stellar wind (~ 1000 km s" 1 ). 

The strong UV resonance lines from Nv, Si IV, and Civ, first seen with IUE (Howarth 
& Wilson 1983b), presumably result from photoionization of the accretion disk and HZ Her 
(Raymond 1993, Ko & Kallman 1994). These lines are much stronger than the strongest 
optical high ionization line, HellA4686, both in absolute flux and equivalent width. Ob- 
servations with the Faint Object Spectrograph (FOS) aboard the Hubble Space Telescope 
(HST) showed that these lines are still present at a few percent of maximum brightness 
during mid-eclipse when the disk and heated star should be entirely obscured (Anderson et 
al. 1994). The source of this emission may be an expanding wind. 

The Goddard High Resolution Spectrograph (GHRS) on HST first resolved these emis- 
sion lines at ps 35 km s _1 resolution (Boroson et al., 1996) to discern variable broad and 
narrow emission components. 

The HST Space Telescope Imaging Spectrograph (STIS) confirmed that the resonance 
lines have at least two components (Vrtilek et al. 2001). A broad component arises on the 
accretion disk while a narrow line component may be associated with HZ Her. Prior to 
the STIS observations, the accretion disk had only been observed as it contributed to the 
continuum light curve or obscured the central X-rays. During eclipse ingress and egress, the 
broad lines seen with STIS behaved as expected for lines from an accretion disk rotating 
prograde with the orbital direction. The blue edge of the line was obscured first in eclipse 
ingress and appeared first in egress. 

Observations with the Far Ultraviolet Spectroscopic Explorer (FUSE) are complemen- 
tary to the existing HST observations. Long-term variability in the system makes it impos- 
sible to combine rigorously data from different epochs. Analyzed separately, however, the 
FUSE wavelength range of 900-1200A offers similar advantages and powerful consistency 
checks to analysis of the HST bandpass of 1200-1700A. Both wavelength ranges have strong 
resonance lines that respond to X-ray photoionization. Observations with the Hopkins Ul- 
traviolet Telescope (HUT) showed that the Ovi doublet has flux comparable to the Nv 
doublet, the brightest near UV line (Boroson et al. 1997). 

The resonance line doublets offer optical depth information through their doublet ra- 
tios, but if the lines are as broad as the doublet separation, it may be impossible to 
determine the individual contribution of each line where they overlap. The separations 
of CivAA1549,1551, SiivAAA1393, 1403, N VAA1238.1242, O VIAA1032, 1038, SviAA933,944 
correspond to Doppler shifts of 500, 1900, 960, 1650, 3600 km s _1 , respectively. The far UV 
lines O VI and S VI compare favorably with the near UV lines (i.e. have greater separation 
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and are less likely to overlap), except for Siiv, which suffers from confusion with an OlV 
blend near 1400A. 

For the present observations, we observed an entire 1.7 day binary orbit with FUSE. A 
major goal of this program was to apply the method of Doppler tomography (§7), which has 
the advantage of diagnosing the accretion flow of different systems without bringing to bear 
more than a few assumptions. 

Table 1 shows the measured physical parameters of the system and parameters we adopt 
for our models. 



2. FUSE Observations 

With FUSE (the Far Ultraviolet Spectroscopic Explorer) we extend observations in the 
UV spectral range to 900-1 190 A, a range observed only once before using HUT, the Hopkins 
Ultraviolet Telescope aboard the ASTRO-1, carried aboard the Space Shuttle but not placed 
in orbit (Boroson et al. 1997). While HUT had a resolution of « 3A, FUSE has a resolution 
of ~ 0.05A. 

FUSE is a NASA Origins mission operated by The Johns Hopkins University. Four 
aligned telescopes feed two identical far-UV spectrographs. With resolution R= 20000 FUSE 
approaches HST in its utility for our program; the time coverage of the HST Space Telescope 
Imaging Spectrograph (STIS) was limited because the detectors were turned off when the 
spacecraft passed through the South Atlantic Anomaly. The FUSE mission is described in 
more detail in Moos et al. (2000) and its on-orbit performance is described in Sahnow et. 
al. (2000). 

Our FUSE observations began on June 9, 2001 at 7:47 UT. Table 2 shows the log of 
exposures, each of which is integrated over each FUSE orbit of the Earth, with gaps when 
Her X-l goes below the horizon. We use the orbital ephemeris of Deeter et al. (1991) to 
determine the orbital phases of our observation. The exposure times listed in Table 2 are 
in most cases equal to the raw observation time. However, for cases such as observation 
13, interrupted by a passage through the South Atlantic Anomaly and not occultation by 
the Earth, the Exptime listed is in general the minimum good exposure time for any FUSE 
detector. The data were obtained through the LWRS aperature and in TIMETAG mode. 

We used the CalFUSE pipeline software version 3.0.8 to extract and calibrate the data 
from all four FUSE telescopes. Below HOOA, where emission features are sharp, we added 
an offset to each wavelength scale, in intervals of 0.025A, so that the absorption and emission 
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features from each detector best agreed. We tested our wavelength calibration against the 
interstellar Si IIA1020.6989 absorption line, which we found to have a mean heliocentric 
velocity of —30 km s^ 1 . The standard deviation of the centroid of this line from orbit to 
orbit was pa 0.03Aor < 10 km s _1 . 

The S/N of the data was pa 5 per 0.1 A pixel in the continuum in the region < 1000A 
and ~ 10 near llOOA. The S/N within the O VI line had greater variation with orbital phase, 
as the doublet changed both in shape and in strength. At = 0.75, the peak S/N within 
the doublet was ~ 15 per O.lA pixel, while at = 0.5, the peak S/N was pa 25 per pixel. 

For further analysis, we skip the region between 1120-1 160 A in the IB LiF detector. 
This region suffers from a systematic decrease in counts known as the worm.. 

In Figures 1 and 2, we show average observed FUSE spectra at orbital phases when the 
emission is dominated by the disk and star, respectively. For the disk-dominated spectrum, 
we use observation number 20, at = 0.917. For the star-dominated spectra, we average 
observations number 7 through 11 (phases = 0.43 0.60). 

For the FUSE spectra at orbital phases dominated by disk emission, we compare with the 
STIS spectrum (observation root name O4V401010) during a Short-On state at = 0.904. 

For the spectrum dominated by star emission, we use the STIS observation with root 
name O4V452050, observed during a Main-On state approximately 5 months earlier, on 
January 24, 2001, starting at MJD 51933.779933 and lasting 2620 seconds. The mean orbital 
phase of this exposure was = 0.455. This spectrum has not previously been published. 

We compare the spectra with the time-weighted average of our continuum models (de- 
scribed in §5). We indicate prominent UV emission lines from the star system. 

Geocoronal (airglow) lines produced in the Earth's atmosphere are also present (see 
Feldman et al, 2001): HlLyman-7 973, Oi5 989, H I Lyman-/3+0 1 4 1026-1027, and Nl 
1167. 



3. Interstellar Lines 

Interstellar features in the FUSE spectra are interesting not only for the direct infor- 
mation they provide on the interstellar medium (ISM), but a proper accounting of these 
features can remove systematic errors from the analysis of the line and continuum emission 
from the system itself. 

The neutral Hydrogen column density, N#, has previously been determined to be ~ 



-7- 



10 20 cm" 2 (Boroson et al. 2000) from the wings of the saturated Lyman a line as observed 
with the HST STIS. The FUSE bandpass includes further saturated lines in the Lyman 
series, and these are consistent with N# ~ 10 20 cm -2 . This N# is also consistent with the 
E(B-V) value of 0.018 according to the Bohlin (1975) relation. 

From the galactic latitude of 37.52° we should expect a sightline with less H and H 2 
than typically seen through the galactic disk. Indeed, while absorption lines from rotational 
levels J = through J = 3 are readily identified, they are not saturated. Thus the level 
populations should be easy to measure and it should be easy to compensate for the effects 
of the absorption lines on the spectra. 

In Figure 3, we show a patch of the time-averaged FUSE spectrum, with the absorption 
profiles expected from source with a flat spectrum absorbed by columns of (2.8, 7.5, 4.1, 
3.4) xlO 14 cm" 2 for absorption from rotational levels J = 0, 1, 2, 3, respectively. We generate 
the H 2 profiles from the templates of McCandliss (2003), which are based on Abgrall et al. 
(1993a,b). We have assumed a line velocity parameter b = 10 km s" 1 and have convolved 
the profiles with a Gaussian to simulate the FUSE Line Spread Function. 

The two lowest rotational energy levels should have populations given by a Boltzmann 
factor, taking into account the statistical weights of the two levels, gi/go = 9. The energy 
difference between the levels is A£ 01 = 170.5 K. We then estimate that the temperature of 
the intervening H 2 gas is T 01 = 140 K. This is hotter than found for H 2 gas from disk stars, 
consistent with the high galactic latitude of the line of sight (Gillmon, Shull, Tumlinson, & 
Danforth 2006). 

The average molecular fraction, defined by 



J N(HI) + 2N(H2) v ' 

is therefore J> 4 x 10~ 5 from the observed H 2 lines, which have N(H 2 )^ 2 x 10 15 cm -2 . Ratios 
of higher rotational levels than J 1 / J are determined not by temperature and collisional 
excitation but by background FUV radiation. The H2 column density observed for Her X-l 
places it just above the boundary at which H2 clouds start to become optically thick to the 
FUV radiation. Gillmon et al. (2006) state that fm ~ 10~ 5 is a typical molecular fraction 
below this boundary. 

Thus the first four levels of rotational excitation of H 2 appear to match observed ab- 
sorption features and to be caused by H 2 clouds that are not out of the ordinary. 
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4. Photometry of Emission Lines 

The lines vary with the binary orbit in a manner similar to that previously observed in 
the near UV with HST. The flux peaks generally near <fi = 0.5, when the X-ray heated face 
of HZ Her points toward the observer, although there may be a dip very close to = 0.5 as 
the accretion disk occults the star. 

In Figure 4 a through i respectively, we show the photometric variation with orbital 
phase of the emission lines S VIA933.4, S VIA944.5, C IIIA977, N IIIA991, O VIA1031.9, O VIA1037.6, 
S IVA1073, P VA1128, and C lllA1176. We have subtracted the background flux using a model 
that we present in §5. For the Nlll line, we remove the flux of nearby airglow lines. 

In Table 3, we show the cross-correlation coefficients between the measured lightcurves 
of the spectral lines and the lightcurves of the disk and star contributions to the continuum 
lightcurve, as determined by our continuum model. 

5. UV Continuum Fits 

The UV continuum varies, as the optical continuum does, with the 1.7 day orbital 
period. The continuum generally peaks near = 0.5, when the X-ray heated face of the 
Roche lobe points towards the viewer, but a portion of the normal star may be blocked by 
the accretion disk, and depending on long-term phase, the actual peak may occur within 
<f>± A0«O.5±O.l. 

The optical/UV orbital light curves do not repeat exactly over the 35-day disk precession 
cycle. The disk may block X-rays from heating portions of the star, and as it precesses it 
projects varying areas into the line of sight. The precessing disk should have different eclipse 
light curves throughout the 35 day cycle. 

To fit the UV continuum observed with FUSE, we assume that X-ray heating of HZ Her 
and the accretion disk cause the entirety of continuum emission. The details of our simu- 
lation are similar to those of Vrtilek et al. (1990, especially the appendix), and Howarth & 
Wilson (1983a) (which was an elaboration on an earlier analysis of Gerend & Boynton, 1976), 
and binary simulation codes by Wilson & Devinney (1971). The method is summarized in 
Appendix A. 

We fit the spectrum using a reddening E(B — V) = 0.018 and for the extinction curve, 
we use Equation 5 in Cardelli, Clayton, & Mathis (1989). 

We also use a more modern determination of the distance to Her X-l (Reynolds et al. 
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1997). This greater distance requires a larger M to reach the same continuum flux. Thus 
the scale of our M values are large compared with those reported earlier from IUE and 
HST observations. The IUE observations found a range of — logM = 8.30 to 8.68 (given 
changes in the long-term phase during an anomalous low state), whereas the HST found 
- log M = 8.19 ±0.06. 

When fit our continuum model to the observed FUSE spectrum from each FUSE obser- 
vation, we find a range of — logM = 8.12 to 8.38. 

For our fits we ignore wavelength regions that contain prominent emission or absorption 
lines. We correct for absorption from the first 4 rotational energy levels of interstellar H 2 by 
multiplying the spectrum by the model presented in §3. 

Table 3 gives the best-fit values of the mass accretion rate, M, which we allow to vary 
as a free parameter for each FUSE orbit. We also list the reduced x 2 - 

We have used our continuum fits to examine the continuum photometry. The model 
allows us to separate the continuum emission into contributions from the heated face of 
HZ Her and from the accretion disk. In Figure 5 we show the model flux near the O VI lines 
split into disk and stellar components. 

6. Eclipse Models 

We made a simple model for emission from a symmetric precessing accretion disk un- 
dergoing eclipse by the Roche lobe of its companion. We then allowed free parameters of 
that model to vary in order to fit the observed O VI profiles during eclipse ingress and egress. 

The disk and Roche lobe geometry are based on the model of Howarth & Wilson (1983a), 
while the formation of the emission lines follows the model of Home (1995), and is described 
in Appendix B. We allow as free parameters the exponent of a radial power law of optical 
depth and elements of the Mach turbulence matrix. We allow the normalization of the line 
flux at each orbital phase to vary as a free parameter. 

We calculate the eclipsing edge of the Roche lobe of HZ Her assuming corotation with 
the orbit. 

Before fitting the model to the data, we subtracted the mid-eclipse spectrum from all 
spectra. The origin of the UV mid-eclipse spectra was explored by Anderson et al. (1994) 
who concluded that it probably does not arise in the accretion disk. The mid-eclipse O VI 
lines are broad, extending from -200 to +500 km s _1 heliocentric velocity. 



- 10 - 



We also apply narrow gaussian absorption lines to our model spectra in order to simulate 
a possible C II interstellar line near 1036A, a possible C II* line near 1037A, and a O I line 
near 1039A. We also include narrow interstellar absorption near the rest wavelengths of 
the Ovi doublet. A narrow absorption line can be seen near the blue component of the 
doublet at A = 1031. 9A, but the weaker line at 1037. 6A cannot be seen explicitly. We fix 
the wavelengths, widths, and optical depths of these lines and do not let them vary in our 
fit. 

We fix all the parameters describing the disk and orbit to those given in Howarth & 
Wilson (1983a), except the outer radius we fix to 2 x 10 11 cm, following Cheng et al. (1995). 

The results are shown in Figure 6. 

The reduced x 2 °f the fit was 2.8 with 1418 degrees of freedom. 

The power law index for the radial dependence on optical depth was a = —1.2, between 
the a = —1.5 result found for the Nv doublet by similar methods in Boroson et al. (2000) 
and the w —0.55 value expected from the simulations of Raymond (1993). 

Although we based our fits on the model of Home (1995) which allows for anisotropic 
turbulence, the fits do not unambiguously settle on particular values of the Mach matrix. 
Here, for simplicity, we present a model with isotropic Mach=l turbulence. 

As with the resonance doublets in the near UV, the observed blue component of Ovi 
near eclipse is stronger than the red component, which suggests that much of the line flux 
is formed where r <^ 1. However, the doublet ratio is uncertain because strong interstellar 
C II and O I absorption lines affect only the red O VI doublet component. Weaker interstellar 
lines, which we have not modeled, may be present as well, and may preferentially absorb one 
or the other O VI component. 

The model lines are double-peaked outside of eclipse, but one peak of the red component 
is absorbed by interstellar O I absorption. The observed line profiles are double-peaked at 
4> = 0.84 and double-peaked but broader at = 0.80. Observations at = 0.35 and = 0.39, 
when narrow line emission predominates, show a dip in the flux at the position of the gap 
between peaks from the accretion disk line. 

The line normalizations we found from our disk eclipse model can be compared with 
the mass accretion rates inferred from our continuum fits. They do not appear correlated. 
At some of the phases investigated, we observe only a small fraction of the accretion disk, 
and none of the illuminated star, so that systematic errors in the line or continuum models 
may make accurate estimates of M impossible. In particular, we have considered the disk 
flat when it is probably warped, we have not calculated the ionization structure of the disk, 
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we have ignored possible emission by the accretion stream, and we have ignored the role of 
winds in this system. 

Chiang (2001) presented models of accretion disk eclipse in Her X-l that considered 
a disk wind. These models were motivated by the lack of double-peaked line profiles in 
any observations of UV resonance lines. Although the current model predicts double-peaked 
emission lines when the full disk is visible, the S/N is not strong enough to rule out a double- 
peaked structure. In addition, interstellar absorption lines coincide with three of the four 
peaks of the O VI doublet. 

It is clear that a model of a partially eclipsed accretion disk with parameters considered 
standard from previous work provides an excellent match to this data set. There appear to 
be too many ways to improve the current fit to justify singling one out. 

6.1. Anomalous emission at = 0.876 

The fit at = 0.876 is particularly poor and we suggest that a process in addition to 
the partial eclipse of a Keplerian accretion disk affects the line profile. 

Figure 6 shows that even though more of the disk is eclipsed than at = 0.835, the 
emission lines at = 0.876 are actually stronger. There also appears to be an additional 
absorption component near the rest wavelength of OviA1032. 

The continuum at = 0.876 has not increased over the continuum at = 0.835 as much 
as the line emission has increased. This is consistent with the earlier result that the X-ray 
illumination from our line and continuum models at other phases were not correlated. 

7. Doppler Tomography 

The Doppler tomography method, developed by Marsh and Home (Marsh & Home, 
1988, Marsh 2005) takes as input an emission line which is broadened by line of sight motion. 
The line must be observed over a good sampling of the binary orbit. In analogy with medical 
tomography, a successful Doppler tomogram builds a higher- dimensional view of the object 
from different observational slices. 

Analysis of Doppler shifted emission lines is prone to the confusion of different gas 
sources which have been projected onto the same observed velocity. The tomographic method 
presents a transform of the same data into features separated not only by projected velocity, 
but also by orbital variability. 
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The output tomogram is not a physical image of an accretion disk, but rather an image in 
"projected velocity space." A component of the emission that varies sinusoidally in velocity 
(but remains constant in flux) is placed at a point in velocity space along a circle with a 
radius equal to the orbital velocity of that emission component. The phase of the sinusoid 
determines its position along that circle. 

The method assumes that all of the variation in a spectral line must be caused by the 
orbit changing the projection of the motion into the line of sight. (For example, the change 
in a spectrum must not be caused by a bright spot being obscured by an opaque object.) 
Even outside of the conditions in which it strictly applies, however, tomographic images can 
provide a common view of the data and may suggest directions for further analysis. 

Technically, a Doppler tomogram is an inversion of 



/oo /*oo 
/ I(v x ,v y )g(v - v R )dv x dv y (2) 
-oo J —oo 

which gives the flux F{y,<f>) at the projected velocity v and the orbital phase as a 
summation of the velocity space tomogram I(v x ,v y ) with a line broadening function g(V — 
Vr) which is usually narrow. This integral is an instance of a Radon transform. The 
line broadening function for our FUSE spectra is approximately 0.05Aand for tomographic 
analysis we assume it is a delta function. 

The relation between the radial velocity and the two axes of the Doppler tomogram, v x 
and v y is given by 

v r — 7 — v x cos 2-7T0 + v y sin27T0 (3) 
where 7 is the systemic velocity of the system. 

The formation of the spectrum from the tomogram can be thought of as follows. The 
phase determines a direction in the velocity space plane, and moving along that line, 
summing the intensity perpendicular to the line, one ideally forms the spectral line. At 
= 0.25, one views from the right in our figures, at = 0.5 one views from the bottom, and 
at = 0.75 one views from the left. 

In this paper, we determine I(v x ,v y ) from F(v,4>) by means of Fourier-Filtered Back 
Projection, which inverts Equation [2] through 

I(v x ,v y ) = / /(7 - cos 27T0 + f y sin 27T0, 0)d0 (4) 
Jo 
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where to calculate f(v , <fi) from f(v, <ft), one takes its Fourier transform, multiplies by a ramp 
filter and a Wiener (optimal) filter based on the noise level, and then takes the inverse Fourier 
transform. 

Although it may seem as if this filtering complicates the method of back-projection, a 
ramp filter (multiplying the Fourier transform by the frequency) is required in order for back 
projection to produce a rigorous inversion of Equation [2] Otherwise, back projection would 
only produce a smeared version of the true tomogram. 

Unfortunately, using a ramp filter also magnifies pixel to pixel noise. We apply a 
standard Wiener filter that assumes a signal amid a noise component that is constant with 
frequency. We find the power spectrum of the line spectra is well fit by several components 
of the form P(u) = Poexp(-kuj) (with k > 0). The contribution of the signal is negligible 
at high frequencies u, where the power spectrum is fit well by a constant noise level n. The 
action of the Wiener filter is then to multiply the power spectrum by P(u)/(n + P(u), which 
is nearly 1 when the signal dominates the noise and negligible at high frequencies where the 
spectrum is almost entirely noise. We used the same form of P(u) for each line at each 
orbital phase, and made a visual comparison between the actual power spectrum and the 
function P(u). Because the filter adapts to different noise levels, the amount of smoothing 
varies from line to line. For the lines with highest S/N, such as OviA1032, the turnover 
frequency for the smoothing filter is on the order of 3 pixels, whereas for weaker lines such 
as S VlA944, the smoothing is on the order of 10 pixels. 

We did not apply the popular method of Maximum Entropy regularization, which awards 
higher probabilities to those tomograms that provide the least information in the sense that 
they are smoothest. 

7.1. Applications of Doppler Tomography 

In Figures 7-14 we show Doppler tomograms of the emission lines of S VI at 933 and 
944A, NlllA911, Ovi at 103lA and 1037A, PvA1128, SivA1073, and C lllA1175. The tomo- 
gram of C IIIA977 is not shown as it is contaminated by the presence of saturated absorption. 
We have interpolated between orbital phases for the integration, which we perform between 
<p = 0.2 and <fi = 0.85. The grayscale of the tomograms extends linearly from the minimum 
to the maximum. 

To understand the Cm multiplet, we consult CHIANTI, a software package and atomic 
database described in Dere et al. 1997 and Landi et al. 2006. The database, used extensively 
by stellar and solar astrophysicists, contains energy levels, wavelengths, radiative transition 
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probabilities, and excitation data for many ions. The associated software package, writ- 
ten in IDL (Interactive Data Language), can be used to examine how line ratios vary with 
temperature and density, subject to certain limiting assumptions. We find that Cm con- 
tains emission components at 0,85,168,199,269, and 367 km s— 1 relative to the 1175. 59A 
component. 

We were able to reduce the artifacts by deconvolving with a triple-gaussian g(V — Vr) 
in Equation [2j We tried Gaussians of equal weight at velocities (0, +168, +367) km s -1 
(Figure 14) and find there is a single bright spot near the Roche lobe. 

A further problem with the Cm multiplet is that the model of the continuum in this 
region predicts an absorption dip. To reduce artifacts, we use the simpler and smoother 
Kurucz model atmospheres, with 1A resolution for all of the tomograms, instead of the 
models based on actual stellar spectra observed with FUSE, which we have computed with 
0.1 A resolution and which include some counting statistics noise. 

The signal in the tomograms, except as noted above, is concentrated in a peak offset 
from the Roche lobe of HZ Her. The traditional signature of an accretion disk is not obviously 
present. Accretion disks are expected to cause broad line emission that appears as a ring in 
a tomogram, with diminished flux within the velocity at the edge of the accretion disk, or 
about 300 km s _1 for Hercules X-l. If the disk emission is symmetric, we expect it to be 
centered on the position of the neutron star in velocity space, (vx,vy) = (0, —169). 

7.2. Narrow Emission and Broad Absorption Components 

There is a simple interpretation of the tomograms in terms of P Cygni lines commonly 
seen in stellar winds from massive stars (and in HST spectra of Her X-l, Boroson, Kallman, 
& Vrtilek 2001). P Cygni lines, caused by resonance scattering, are characterized by red- 
shifted emission and absorption blue-shifted by velocities common in the wind. There is a 
hint of blue-shifted absorption in the trailed spectrograms at = 0.3 — 0.6, although there 
may be an interstellar absorption line as well at these wavelengths and the gap between the 
peaks of the accretion disk spectrum could also appear to be absorption. 

We note that for many tomograms in addition to the bright spot at (vx, vy) ~ (150, —30), 
there is a dark spot at (vx,vy) ~ (—100,50). (For OviA1031, this absorption may have a 
darker lobe near the rest wavelength from interstellar Ovi absorption.) Near the start of 
the phase range used by our tomograms, at <fi ~ 0.3, the observer's line of sight comes from 
the lower right in Figure 7. As the line of sight passes through both bright and dark spots, 
the observer at <fi ~ 0.3 sees a diminished central peak in the spectral lines. 
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At ps 0.5 — 0.6, the viewer looks at the system from the bottom of the plot, stack- 
ing pixels vertically. The result is red-shifted emission (the bright spot) and blue-shifted 
absorption (the dark spot). 

Although the blue-shifted absorption at = 0.3 — 0.6 may be explained through the 
analogy of P Cygni lines, the red-shifted emission at these phases does not fit this explanation. 
In a P Cygni line, except when a very hot wind is present (as in a Wolf-Rayet star), all 
the emission is caused by spherically symmetric scattering, so that conservation of photons 
requires that the forward-scattered flux be nearly equal to the back-scattered ( "absorbed" ) 
flux (some portion of the forward-scattered emission is blocked by the star itself). In contrast, 
any blue-shifted absorption in the Her X-l spectrum is only a small fraction of the emission. 
Further, P Cygni emission is only red-shifted because it overlaps with the absorption on the 
blue end; the actual emission spectrum should be symmetric about zero velocity. 

It seems most natural to associate the bright spot with the X-ray heated face of the 
normal star, as it follows roughly the same light curve and appears near the star and far 
from the disk on the tomograms. Yet the placement of the spot, or equivalently, the orbital 
variation with velocity, cannot be accounted for easily by such a model. 

In Figure 15, we show the observed spectrum near Ovi (bold), the spectrum expected 
if the illuminated surface of HZ Her emits Ovi Doppler shifted by the local rotational 
velocity (dashed), and an empirical Gaussian model of the emission lines. The Doppler 
shifts predicted from the rotation of the surface of HZ Her are less than observed in the 
narrow emission lines. 

To model the narrow O VI emission lines we adapt our model of the continuum emission. 
This model takes into account the X-ray shadow cast by the disk and the eclipse of portions 
of the normal star by the disk. The reprocessed O VI emission is proportional to the incident 
X-ray flux. We also experimented with a model in which the X-ray illumination was non- 
isotropic, illuminating HZ Her with greater luminosity at = 0.25 than at = 0.75. This 
model still does not account for the observed velocities. 

We also made an empirical model of the lines as Gaussians. The flux in the Gaussian 
line emission at each phase was proportional to the flux in the stellar component of our 
continuum model using Kurucz model spectra. The central velocity in this model varies 
sinusoidally with phase. We added Gaussian absorption to correspond to features seen in 
HST STIS spectra of the N VA1240, SilvA1400, and CivA1550 doublets near = 0.3 - 0.6. 
We set the maximum projected redshift of the emission lines and blueshift of the absorption 
to occur near = 0.5, but delayed by = 0.06, corresponding to a deflection of a wind 
by the Coriolis effect. The absorption line has maximum covering fraction near = 0.5 
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that smoothly varies to at <fi = 0. (Any wind is probably confined to the cylinder from 
the star towards the disk, as the emission lines seen in mid-eclipse are only ~ 1% of the 
peak line strength, as observed with the HST FOS by Anderson et al. 1994). Both emission 
and absorption components have constant width. Both emission and absorption components 
vary about a central velocity given by the LI point. 

When these features are added to the model of the disk spectrum in §6, assumed to 
move with the neutron star's 169 km s _1 orbit, and random counting statistics noise is added, 
the resulting tomogram (Figure 16) resembles the observed tomogram, while showing only 
a hint of the presence of an accretion disk. Without the assumed phase delay of A0 = 0.06, 
the bright spot would appear rotated to vy — 0. 

8. Line Ratio Diagnostics 

The FUSE spectra show emission lines that may serve as diagnostics of density, tem- 
perature, or optical depth. 

8.1. C ill line diagnostics 

The simultaneous measurement of ClllA977 and ClllA1176 can provide us with an 
important diagnostic. The C ill 1176/977 ratio is sensitive to density in the range up to 
about n e = 10 11 cm" 3 . At higher densities (up to n e ~ 10 15 ), the populations of the ground 
and metastable levels are entirely determined by collisions. In that regime, the 1176/977 
ratio depends only on temperature and the optical depths in the lines. 

We ignored flux from the FUSE LiF2B detector near the C IIIA977 line, even though 
this detector recorded flux down to exactly 977A. The edge of the detector causes the flux 
measured in this region was significantly lower than the other two detectors, so that including 
this data would have created an artificial dip at wavelengths A > 977A. 

The fluxes we measure for C IIIA977 are lower limits because of a saturated interstellar 
absorption line at the rest wavelength. We correct the fluxes for interstellar reddening using 
E(B-V)=0.018. 

We compare the fluxes of the two lines near = 0.6, when the fluxes peak. If the 
C IIIA977 line behaves similarly to the O VI lines, then the bright narrow emission will, near 
this phase, have the greatest redshift. If so, the line may be redshifted away from the 
interstellar absorption and the flux measurement may be more accurate near this phase. If 
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the O VIA1032 line were subject to similar interstellar absorption lines as the C IIIA977 line, 
the flux at = 0.6 would be diminished by ~20%. 

If we measure the raw ClllA977 flux from the two phases surrounding = 0.6 to be 
7 ± 1 x 10~ 13 erg s" 1 cm" 2 , then we find a 1175/977 ratio of 0.82 ± 0.16. If we measure 
the ClllA977 flux at exactly the peak at = 0.6 and attempt to correct for interstellar 
absorption lines, we find a ratio of 0.60±0.09. 

These ratios are at the high end of those predicted by the CHIANTI software, which uses 
atomic data from the CHIANTI database and computes line ratios, given certain simplifying 
assumptions (for example, the lines are optically thin and the gas is collisionally ionized and 
not photoionized). If we consider a density of n e = 10 13 ' 4 cm -3 , the lower ratio found above 
would restrict the temperature to T > 6 x 10 4 K. 

The ratio is rendered uncertain observationally by the saturated interstellar absorption 
line, and theoretically because of optical depth effects. Raymond (1993) presented models of 
line emission from an X-ray illuminated accretion disk which took into account optical depth 
through an escape probability formalism. One version of the model, "COS", assumed cosmic 
abundances while the other, "CNO", assumed that abundances had been altered by CNO 
processing. The CNO models should be more appropriate for Her X-l. In the two cases the 
J(1176)//(977) ratios were 1.3 and 1.1, respectively. The observed emission at = 0.6 may, 
however, have an origin not on the disk but on the normal star. 

8.2. Bowen fluorescence and NlllA991 emission 

The Bowen Fluorescence process (Schachter et al., 1989) arises because of the nearly 
perfect coincidence of the Hell Lya and the O ill 2p 2 - 2p3d resonance line (A304), resulting 
in Olll near-UV primary cascades at AA3133, 3444 and secondary cascades at A374 back 
to the ground state. If conditions are right, an additional fluorescence occurs, since the 
O IIIA374 line is almost coincident with the two N III 2p - 3d resonance lines, resulting in 
N III optical primary cascades at A4634, 4641, 4642. The emission in all of these cascades is 
completely dominated by the Bowen process; detection of any of them is a clear confirmation. 
We note that in Her X-l the process requires substantial optical depth (~ 10 6 ) in the 
He II Lya and O III A374 pumping lines, because otherwise they will simply escape without 
conversion to a Bowen line. 

In the Sun, Raymond (1978) discovered O IIIA304 Bowen emission. It was found that 
measurements of O IIIA703, a ground-state Bowen cascade which, therefore, can also be 
produced by collisional excitation, can serve as a density diagnostic, and compared favorably 
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with other solar estimates of n e . 

We may use NlllA990, an analogous ground-state Bowen cascade produced as a result 
of the N III primary cascades ( "A4640" ; all other III and N III ground-state cascades lie 
below the Lyman limit and hence are unobservable). 

The intensity of N IIIA990 resulting from Bowen fluorescence is related to the intensity 
of the Bowen lines at 4640Aby 

7/(991) = (4640/991)5(3p2p)7(4640) (5) 

where 7?(3p2p) is the branching ratio of 2p-2s 2 p 2 2D versus 3p-3s. From CHIANTI, we take 
fi(3p2p) = 0.55 and 5(3d3p) = 0.0058. 

The remainder of the line intensity is therefore from collisional excitation and is given 
by 7 C (991) = 7(991) - 7/(991). 

Taking into account the reddening toward Her X-l implied by E(B-V)=0.018, we find a 
peak flux in the 991A line of 8.6±0.7x 10~ 13 erg s -1 cm -2 . Non-simultaneous measurements 
of the optical A4640 line have a peak flux at similar orbital phases of 7.0 ± 0.2 x 10~ 14 
erg s _1 cm" 2 (Still et al. 1997). The optical flux varies from orbit to orbit, and, moreover, 
probably includes contributions from C ill as well as N ill. 

These values imply that 20% of the N IIIA991 line flux is the result of the Bowen fluo- 
rescence mechanism. 

We attempt to relate the density to the Bowen flux by 

4 (6) 

n e N(Nm)q 
B(3d3p)5(3p2p)iV(N iii)o-7(374) 
4.9 x 10 8 n e /7(374) 

with iV(Niii) the density of Nm and q = 6 x 10~ 8 s _1 the excitation rate of the 99lA 
transition at T=30,000 K. 

We estimate an upper bound on the O IIIA374 intensity (the line may be optically thick) 
from the O IIIA3132 line observed with the Faint Object Spectrograph on HST by Anderson 
et al. (1994). We account for the other Om branches leading to the A374 line by using the 
observed Om Bowen spectrum of RR Tel (Selvelli, Danziger, & Bonifacio 2007). We then 
estimate 7(374) <^ 0.032(7J)/_R) 2 , with D = 6.6 kpc the distance to the system and R the 
radius of the emitting region. Boroson et al. (2000) found logn e = 13.4 ± 0.2 for the narrow 



7 C (991) 
77(991) 
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lines from the ratios of lines observed with HST, and Howarth & Wilson (1983b) found 
\ogn e = 13.3 from IUE observations. Combining with Equation |6j we find R < 7 x 10 10 cm. 

This is smaller than the size of the accretion disk (-Router ~2x 10 11 cm). We note that 
the tomograms show emission from a small region in velocity space. 

Another test is provided by a comparison with theoretical models of NlllA991 emis- 
sion in the absence of the Bowen process. Raymond (1993) found ratios between the Nlll 
line and ClllA977 for the COS model of I(991)/I(977) = 0.34 and for the CNO model of 
J(991)/J(977) = 0.53. The observed ratio of 0.75 ± 0.10 requires enhancement by Bowen 
fluorescence and provides an excellent match to the CNO model if only 70% of the N IIIA991 
line results from collisional excitation. 



9. Discussion 

The far UV spectrum of Hercules X-l shows clear evidence for a Keplerian accretion 
disk. A simple model such a disk can fit the broad Ovi emission lines near eclipse. The 
origin of brighter, narrow emission lines is still unknown. Doppler tomograms place this 
emission apart from the Roche lobe of HZ Her. 

The bright narrow emission component generally follows the flux expected from the 
illuminated portion of the normal star. However, it appears brighter than expected at 
= 0.2 and the velocity excursion is also greater than is accounted for by the rotation of 
the Roche lobe. A dense portion of a stellar wind that moves along the surface of the star, 
or that originates at the LI point and flows back towards the star, would be red-shifted 
at cf) = 0.5 and blue-shifted at = 0.25 and <fi = 0.75. Although such a model can be 
made to reproduce the gross features of the line variation, there remain too many possible 
alternate explanations for this to be convincing. There could be a source of variable emission, 
absorption, or scattering in addition to the normal star, for example, the gas stream or the 
surface of the disk. 

In the presence of stellar winds, possibly transient or restricted in solid angle, tomogra- 
phy did not show a clear signature of an accretion disk in Hercules X-l. 

Although we found that a simple model of an accretion disk generally fit the spectra 
during eclipse ingress and egress, there was an anomalous observation in which the O VI lines 
brightened as the disk eclipse progressed. The phase <fr — 0.876 is in the <fi = 0.8 — 0.9 range 
in which historically "pre-eclipse dips" of X-rays have been observed, and it is tempting to 
relate the temporary brightening of the lines to a physical cause or result of the dips. 
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Alternately, an anomalous change to the emission lines could result from a change in 
shadowing of the X-ray emission. Boroson et al. (2001) concluded that anomalous absorption 
within the N VAA1238, 1242 doublet changed too rapidly to result from material passing over 
the line of sight. The absorption changed so rapidly probably because of a rapid change in 
the shadowing of the X-rays which ionize the gas. 

Future X-ray spectroscopy missions such as Constellation X may be able to resolve lines 
from more highly ionized species than the UV and Far UV resonance lines (Vrtilek et al. 
2004). If more highly ionized species track disk material better and are not as prevalent in 
the narrow line region, we may be able to make reliable tomograms of the disk in this or 
similar systems. Hydrodynamic models of this system may be required to cut down on the 
vast parameter space of possible gas flows responsible for the spectral signatures observed 
here. 

We would like to thank Jeff Bryant for helping us develop some Mathematica routines for 
Doppler tomography. CHIANTI is a collaborative project involving the NRL (USA), RAL 
(UK), MSSL (UK), the Universities of Florence (Italy) and Cambridge (UK), and George 
Mason University (USA). We would like to thank the anonymous referee for suggestions 
that improved the text and presentation, for pointing us to the SviAA933,944 lines, and 
prompting us to analyze the interstellar molecular absorption lines. 

A. Details of the Model of the Her X-l Continuum 

With one free parameter, the mass accretion rate M, we model the continuum emission 
from the accretion disk and X-ray illuminated face of the normal star. We do this not 
by means of radiative transfer calculations, but by calculating the heated temperatures of 
both surfaces, and co-adding spectra appropriate for those temperatures, either blackbody 
spectra, models of stellar spectra, or actual spectra of hot stars observed with FUSE. 

Our model of the far UV continuum emission of Her X-l follows closely the methods of 
Vrtilek et al. (1990) and Cheng, Vrtilek, & Raymond (1995). We describe the method here 
in detail, including the minor departures we have made to extend the spectral simulation to 
the far UV. 

The model of the accretion disk temperature as a function of radius, T(r), includes 
the effects of both internal heating due to viscous forces and heating from X-rays from the 
neutron star. 

When we calculate the disk temperature as a function of radius, T(r), we include both 
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the heat from viscous forces and X-rays from the neutron star. The local energy generated 
by accretion is given by 

4 3GMM I IR 

aT o( r ) = _„ 3 1 " \ T ( Al ) 



87rr 3 \ V r 



for Stefan-Boltzmann constant a, neutron star mass M and radius R, and disk mass accretion 
rate M (Shakura & Sunyaev 1973). 

For local energy balance, the energy emitted must equal the energy generated plus the 
X-ray energy absorbed: 

"iw^*)^^ (A2) 

where we choose the albedo a = 0.5, and the disk height at radius r is h(r). The X-ray 
luminosity is related to the gravitational potential released by accreting matter: 

GMM , K x 

L x = 0.5—— (A3) 
a 

These equations, together with vertical hydrostatic equilibrium, can be solved numer- 
ically to determine T(r). The surface of the noncollapsed star is also assumed to have an 
albedo a = 0.5. 

At each small region in the accretion disk, the temperature T(r) is compared with the 
critical disk temperature, T dc . If T(r) > T dC) then we assume that the disk region radiates 
like a blackbody. If T{r) < Td c , then we interpolate through our library of actual stellar 
spectra to find a spectrum appropriate for temperature T{r). From experience fitting HST 
FOS spectra and fitting the Balmer jump, we fix T dc = 10000 K so that in practice the 
disk almost always has blackbody spectra, except at its edge where it is not illuminated by 
X-rays. 

The stellar library from 1150 to 7500A is described in Cannizzo & Kenyon (1987). The 
near UV spectra were obtained with IUE and are described in Wu et al. (1982). All the 
spectra are normalized to have constant flux in the V bandpass and are then scaled to their 
absolute magnitudes using their V-R colors and the Barnes-Evans relation (Barnes, Evans, 
& Moffett 1978). 

For this paper, we have extended the library of stellar spectra further into the UV by 
using FUSE spectra of main sequence stars between 9.5 and A 7. These observations are 
described briefly in Table 5. Stars of temperature greater than T sc = 18, 900 K are not used 
in the fits to Her X-l and are not included in the table. We have interpolated over those 
wavelengths we expect are affected by interstellar atomic or molecular absorption, and we 
have applied a filter designed to smooth interstellar lines while retaining stellar features. We 
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de-reddened each spectrum and scaled the flux according to the star's visual magnitude V. 
We have tested the method and the use of these stellar spectra by using model stellar spectra 
as predicted by Kurucz (1979). Plots of far UV spectra in our library, normalized to have a 
constant V magnitude, are shown in Figure 17. 

For the purpose of computing the disk shadow, the disk is assumed to have a fixed 
opening angle and a fixed tilt from the orbital inclination, although the direction of the disk 
normal precesses with the X-ray cycle in a manner described by Gerend & Boynton (1976) 
and Howarth & Wilson (1983a). Because of the finite disk opening angle, the outer regions 
of the disk can occult the inner regions, and the disk edge itself can radiate, although it is 
not illuminated by the central X-rays. 

The noncollapsed star is assumed to fill its Roche lobe, and its shape, visibility, and 
illumination by the X-ray source are treated according to methods evolved from Wilson & 
Devinney (1971) to Howarth & Wilson (1983a). Following those references, we model the 
eclipse of the disk using the Roche potential. We take into account both limb darkening and 
the gravity-darkening appropriate for a late-type star. If the temperature of some region on 
the surface of the star, after X-ray illumination, is greater than some critical temperature T sc , 
we assume the region emits as a blackbody. If the local temperature is less than T sc , we again 
interpolate through our stellar library to find the appropriate stellar spectrum with which 
the spot radiates. From earlier experience fitting Her X-l UV spectra, we fix T sc = 18900 
(Cheng, Vrtilek, & Raymond 1995). 

The model, as currently realized, does not account for emission or shadowing by a gas 
stream between the noncollapsed star and disk. 

B. Details of the Model of the Disk Lines in Eclipse 

Accretion disks in theory have emission lines with a double-peaked shape, with peaks 
separated by a velocity given approximately by the projection into the line of sight of the 
orbital velocity at the edge of the accretion disk. For Hercules X-l, with typical neutron 
star mass, an nearly edge-on inclination, and an outer disk radius of R ~ 2 x 10 11 cm, we 
expect the peaks to occur at ±300 km s _1 . While greater velocities are found within the 
outer disk radius, each ring at a fixed radius emits down to a Doppler velocity of where 
the gas moves tangential to the line of sight. Regions of high Doppler velocity arise from 
progressively smaller regions in the disk, and contribute less to the overall line shape. 

We have calculated detailed models of emission from an accretion disk and have fit these 
to the profiles observed in Her X-l during eclipse ingress and egress. In reality, the disk is 
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likely to be warped, and in future work we will test individual models of this warped shape. 
The importance of the demonstration in this paper is that it shows how easy it is to get 
general agreement with the observed line profiles given the standard picture of the accretion 
disk. 

The method in detail follows that of Home (1995), which allows for nonisotropic turbu- 
lence and for Keplerian shear, that is, the local dispersion in disk velocity as a result of the 
Keplerian flow itself. This introduces free parameters for the "Mach matrix" describing the 
turbulent flow, in addition to the power law exponent that we use to as a phenomenological 
description the strength of the emission line as a function of radius in the disk. For simplicity, 
we keep the elements of the Mach matrix constant with radius in the disk. 

Following Home (1995), the local line profile is given by 

I v oc (1 - e- T ") (Bl) 

where the optical depth r„ at frequency v is given by 

r„ = roe (^o) 2 /2Av- (B2) 

for a given line center optical depth To- 

The nonisotropic turbulence and shear enter through AV, given by 

AV 2 = AVLrm + AKheax + AV** (B3) 

or 

AV 2 /C 2 S = — A + (Q sin i tan i sin 26) 2 + sin 2 i(M RR cos 2 9 -2M m cos 9 sin 9 (B4) 
•~yA 



+ M ee sin 2 9) + 2 sin i cos i(M ZR cos 9 - M ze sm9) + cos 2 iM zz 



(B5) 



where i is the orbital inclination, /i is the mean molecular weight, 7 = 5/3, and A is the 
atomic weight. The shear parameter Q is given by 

3 H Vkep AZ 

Q ~4R^72H (B6) 
and, for simplicity, is set to 3/4, the maximum value for an optically thin accretion disk. 

The Mach matrix is given by the correlations between components of turbulence in 
different cylindrical coordinates: 

M„ S » (B7) 



S 



where C s is the sound speed. 
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Symbol 


Adopted value 


Meaning 


Reference 


Q 


0.58 


Binary mass ratio 


Howarth & Wilson 1983a 




1.4 


Mass of neutron star (M ) 




dinner 


10 6 


Inner radius of accretion disk (cm) 


Cheng, Vrtilek, & Raymond 1995 


-Router 


2xlO n 


Outer radius of accretion disk (cm) 


CVR 


a 


0.5 


Albedo 


CVR 


D 


6600 


Distance to system (pc) 


Reynolds et al. 1997 


Eb-v 


0.018 


Reddening parameter 


Boroson et al. 2001 


T* 


8100 


Polar temperature of normal star (K) 


CVR 




0.8105 


Disk precession phase offset 


H&W 


d 


4.87 


Disk thickness (degrees) 


H&W 


a d 


28.72 


Disk tilt from orbital plane (degrees) 


H&W 


a* 


6.35xlO n 


Orbital separation, centers of mass (cm) 


H&W 


i 


81.25 


Orbital inclination (degrees) 


H&W 



Table 1: Parameters describing the X-ray binary Hercules X-l 
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Root name 


MJD (start) 


MJD (end) 


Orbital Phase 


Exptime (s) 


B0080101001 


52069 


.32464 


52069. 


.36421 





.182562 


3418 


B0080101002 


52069. 


,39311 


52069. 


43359 





.223105 


3498 


B0080101003 


52069 


,46245 


52069. 


.50297 





r\ s~i c\ r\ r\ r\ 

.263900 


3501 


B0080101004 


52069 


,53186 


52069. 


.57233 





.304712 


3497 


B0080101005 


52069 


,60135 


52069. 


.64170 





O A V V A f~7 

.345547 


3455 


B0080101006 


52069 


.67566 


52069. 


71106 





.387799 


2800 


B0080101007 


52069. 


.74947 


52069. 


.78043 





a c\r\r\r\ a 

.429904 


2675 


B0080101008 


52069 


.82233 


52069. 


.84980 





A 171 "TO A 


2374 


B0080101009 


52069 


OA /I An 

.89492 


52069. 


.91917 





.513482 


2095 


B0080101010 


52069 


.96775 


52069. 


.98854 





.555305 


1796 


B0080101011 


52070. 


.04000 


52070. 


.05792 





.596954 


1548 


B0080101012 


52070. 


.11058 


52070. 


.12730 





.638113 


1442 


B0080101013 


52070. 


,15631 


52070. 


.16500 





.662651 


644 


B0080101014 


52070. 


.18235 


52070. 


19666 





.679620 


1231 


B0080101015 


52070. 


,22613 


52070. 


26604 





.712899 


3447 


B0080101016 


52070. 


.29521 


52070. 


.33540 





.753612 


3472 


B0080101017 


52070 


.36448 


52070. 


.40475 





.794380 


3479 


B0080101018 


52070. 


.43388 


52070. 


.47413 





.835192 


3478 


B0080101019 


52070. 


.50334 


52070. 


.54349 





.876021 


3469 


B0080101020 


52070. 


.57267 


52070. 


,61286 





.916809 


3434 


B0080101021 


52070. 


,64515 


52070. 


68223 





.958526 


3204 


B0080101022 


52070. 


,71889 


52070. 


.75160 





.000611 


2811 


B0080101023 


52070. 


.79200 


52070. 


.82097 





.042515 


2503 


B0080101024 


52070. 


.86463 


52070, 


.88778 





.083521 


2000 



Comments 



Dip in VI (double-peak gap 
Dip in VI 



OviA1032 double-peaked 
OviA1032 double-peaked 
Lines anomalously bright 



Mid-Eclipse 



Table 2: A log of the FUSE Her X-l observations. 
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Exposure Orbital Phase Duration (s) M (Kurucz) X 2 (Kurucz) M {FUSE) xl {FUSE) 









(- lo Sio M o) 




{- lo Sio M o) 




1 


(J.loz 


O A 1 O 

o41o 


o n a 

o.(J4 


O 7 

z.7 


o nr 

8. (Jo 


O 7 

z.7 


2 


n ooo 


O A C\0 

o49o 


o 

o.(J7 


o n 

z.y 


o i n 

8. ID 


o n 
O.L) 


6 


n 0£? /i 
U.Z04 


o n 1 
OOUl 


o nn 

o.uy 


O 1 

O.l 


o i n 

o.iy 


o o 
6.2 


A 

4 


n ontx 
U.oUo 




O 1 P 

o.lo 


o n 

z.y 


O OP 

O.ZO 


o n 
o.U 


r 




n o a a 
U.o4o 


O A OP 

o4oo 


o oi 

o.zl 


O 7 

Z. f 


O OO 

o.oz 


o n 

z.y 


p 



n ooo 

0.388 


on rr n 

3059 


O O A 

o.z4 


o c 

z.o 


o o o 

8.33 


O 7 

z.7 


7 


n a on 

U.4o(J 


zo7o 


o on 

o.z9 


O /I 

z.4 


O O O 

o.oo 


o o 
Z.O 


o 
O 


n /i to 


ZO (4: 


o on 
o.zl) 


o o 
z.o 


o on 
o.oU 


O E 

Z.O 


n 

y 


n ci o 
U.Olo 


onn e 

zUyo 


o ri7 
o.U/ 


O 

z.o 


O 1 o 

o.lo 


O 7 

z. f 


1 n 

ID 


(J. 555 


1796 


o no 

8.08 


O A 

2A 


O 1 7 

8.17 


O 7 

z.7 


ii 


(J. 597 


i r /in 

1548 


O 1 o 

8.13 


o c 

2.5 


O O 1 

o.zl 


O P 

2.6 


12 


n poo 

0.638 


1444 


o on 

8.20 


2.2 


O 07 

8.27 


o o 

2.3 


13 


0.663 


751 


8.22 


1.2 


8.30 


1.3 


14 


0.680 


1236 


8.19 


1.6 


8.27 


1.6 


15 


0.713 


3448 


8.17 


2.5 


8.27 


2.6 


16 


0.754 


3473 


8.20 


2.5 


8.22 


2.6 


17 


0.794 


3480 


8.18 


2.3 


8.19 


2.3 


18 


0.835 


3478 


8.22 


2.2 


8.22 


2.2 


19 


0.876 


3469 


8.17 


2.2 


8.18 


2.2 


20 


0.917 


3472 


8.29 


1.7 


8.28 


1.7 


21 


0.958 


3204 


NA 


NA 


NA 


NA 


22 


0.001 


2826 


NA 


NA 


NA 


NA 


23 


0.042 


2503 


NA 


NA 


NA 


NA 


24 


0.084 


2000 


8.10 


1.7 


8.12 


1.7 



Table 3: Fits to the Her X-l continuum as observed with FUSE. M (Kurucz) and \ 2 (Kurucz) 
give the mass accretion rate and goodness of fit parameter for the model fit that uses Kurucz 
spectra while M (FUSE) and x 2 (FUSE) give those parameters for fits using a library of 
FUSE spectra. 
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Ion 


Wavelength (A) 


'"star 


rdisk 


'"continuum 


Comments 


Svi 


933.4 


0.88 


0.38 


0.67 




Svi 


944.5 


0.88 


0.54 


0.85 




Cm 


977 


0.90 


0.45 


0.76 


Affected by strong ISM line 


Nm 


992 


0.80 


0.62 


0.91 


Near strong airglow lines 


Ovi 


1031.9 


0.74 


0.75 


0.95 




Ovi 


1037.6 


0.82 


0.71 


0.98 




Siv 


1073 


0.75 


0.63 


0.88 




Pv 


1128 


0.75 


0.58 


0.86 




Cm 


1175 


0.82 


0.65 


0.95 





Table 4: Correlation coefficients between the light curves of emission lines and continuum 
components as separated by our model (§5). 



Star 


Spectral Type 


T eff (K) 


FUSE Rootname 


E(B-V) 


V 


HD 146813 


B1.5 V 


24000 


P1014901 


0.02 


9.06 


HD 74662 


B3 V 


18000 


A1290201 


0.09 


8.87 


(Interpolated) 


B4 V 


16500 








HD 92288 


B6 V 


14000 


Z9012801 


0.05 


7.9 


(HD 21672,HD 92536) 


B8 V 


11500 


(Z9011401,Z9012901) 


(0.09,0.04) 


(6.63,6.32) 


HD 149630 


B9 V 


10800 


B0910101 


0.04 


4.2 


HD 109573,HD 181296) 


AO V 


10000 


(B0910401,P2500101) 


(0.0,0.01) 


(5.78,5.03) 


HD 31647 


Al V 


9300 


C0380901 


0.0 


4.989 


HD 115892 


A2 V 


9050 


A0410505 


0.0 


2.75 


HD 43940 


A3 V 


8850 


C0380101 


0.0 


5.87 


HD 11636 


A5 V 


8500 


A0410101 


0.0 


2.64 


(Interpolated) 


A6 V 


8350 








HD 187642 


A7 V 


8050 


D0990101 


0.0 


0.77 



Table 5: FUSE spectra used to form a library of stellar continua versus effective temperature 
in the far UV. 
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svisvi c 111 isi 111 ovi sivsiv pv cm nv 




900 950 1000 1050 1100 1100 1150 1200 1250 1300 

Wavelength (A) Wavelength (A) 



OV Si IV NIV CIV Hell 




Fig. 1. — The mean Her X-l far UV spectrum observed with FUSE (< 1180A) and the 
HST STIS (> 1150A) near orbital phases when the emission is expected to be dominated 
by the accretion disk. Prominent emission lines are indicated by dotted vertical lines. We 
show he average of our continuum models using actual FUSE spectra of hot stars together 
with blackbodies. The component of the model resulting from the X-ray heated face of the 
normal star is indicated by a dotted curve. 
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OV Si IV NIV CIV Hell 




Fig. 2. — The mean Her X-l far UV spectrum observed with FUSE (< 1180A) and the HST 
STIS (> 1150A) near orbital phases when the emission is expected to be dominated by the 
X-ray heated face of the normal star. The spectrum is plotted as in Figure 1. The lines and 
continuum are brighter at these phases, and the lines are narrower. 



- 33 - 



Flux (10 13 erg s 1 crrr 2 A -1 ) 



CD O 
00 
O 
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Fig. 3. — A sample of the H 2 absorption lines in the spectrum of Her X-l, and a model 
(dotted lines) with flat normalization, of absorption from rotational levels J = through 
J = 3. 
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Fig. 4. — Light curves of prominent lines in the FUSE spectrum of Hercules X-l. A model 
of the continuum flux has been subtracted. Errors at the la level are shown. Emission 
lines shown: (a) SviA933.4, (b) SviA944.5, (c) CinA977, (d) NlllA991, (e) OviA1031.9, (f) 
OviA1037.6, (g) SivA1073, (h) PvA1128, (i) CiiiA1176 
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Star continuum. 



Total 



a 0.2 a. a a.6 0-S 
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Fig. 5. — Flux versus orbital phase in our model of the continuum flux in the wavelength 
range of the Ovi doublet. The continuum is separated into disk and stellar components. 




Fig. 6. — Fits to the Ovi lines at phases when the accretion disk is partially eclipsed. 
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Fig. 7. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
S VIA933 from = 0.2 to <fi = 0.85. The grayscale runs from the minimum to the maximum. 
Top: the trailed spectrogram (left) and filtered trailed spectrogram (right). Middle: the back- 
projected tomogram (left) and filtered tomogram (right). Bottom: the trailed spectrogram 
inverted from the back-projected tomogram (left) and from the filtered tomogram (right). 




Fig. 8. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
S VIA944 presented as in Figure 7. 
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Fig. 9. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
NinA991.6 (Bowen process) presented as in Figure 7. 
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Fig. 10. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
OviA1031.9 presented as in Figure 7. 
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Fig. 11. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
OviA1037.6 presented as in Figure 7. 




Fig. 12. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
SivA1073 presented as in Figure 7. 
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Fig. 13. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
PvA1128 presented as in Figure 7. 
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Fig. 14. — A Fourier-filtered back-projected Doppler tomogram using the FUSE spectra of 
ClllA1176 presented as in Figure 7. 
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Fig. 15. — The OviA1032 emission line averaged over every two FUSE exposures (thick 
curve), a model of Ovi emission on the surface of the noncollapsed star (dashed), and an 
empirical model of Gaussian emission and absorption models. Both models include flux from 
a model of the accretion disk. 
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Fig. 16. — A Doppler tomogram of a simple simulation of the OviA1032 line as a sum of 
moving Gaussians and an accretion disk model, with counting statistics included. 
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Fig. 17. — Sample spectra, with spectral type labelled, in the library of FUSE spectra we 
use to complement our library of WE and optical spectra. We use this library in the fitting 
of the Her X-l continuum. All stars have their V band flux normalized to the same value. 



